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We show that both DC and AC Josephson effects with superfiuid Fermi atoms in the BCS-BEC 
crossover can be described at zero temperature by a nonlinear Schrodinger equation (NLSE). By 
comparing our NLSE with mean-field extended BCS calculations, we find that the NLSE is reliable 
in the BEC side of the crossover up to the unitarity limit. The NLSE can be used for weakly- 
linked atomic superfluids also in the BCS side of the crossover by taking the tunneling energy as a 
phenomenological parameter. 



I. INTRODUCTION 

Few years ago the AC Josephson effect 0, Q with atomic superfluids was predicted Q and observed Q with 
Bose-Einstein Condensates (BECs). AC Josephson oscillations in superfiuid atomic Fermi gases have been considered 
theoretically by several authors [E IE 0> HI • Recently, Spuntarelli, Pieri and Strinati @ have studied the DC Josephson 
effect across the BCS-BEC crossover with neutral fermions by using the extended BCS equations: they computed 
the current-phase relation throughout the BCS-BEC crossover at zero temperature for a two-spin component Fermi 
gas in the presence of a barrier [9( . 

In this paper we discuss a nonlinear Schrodinger equation (NLSE) [E [Tol . [Tl| which is equivalent to the equations 
of superfiuid hydrodynamics [l2| with the inclusion of a gradient term of the von Weizsacker type [l3L IbH . [l5| . By 
using the NLSE we reproduce the DC Josephson results of Spuntarelli, Pieri and Strinati 0, in the BEC side of 
the BCS-BEC crossover, i.e. from the unitarity region to the deep BEC regime. However, we show that our NLSE 
describes Q both DC and AC Josephson effects [Ell] also in the BCS regime of two weakly-linked atomic superfluids 
if the tunneling energy is taken as a phenomenological parameter a la Feynmann [iH \vl\ . 

II. HYDRODYNAMICS OF FERMI SUPERFLUIDS AT ZERO-TEMPERATURE 

At zero temperature the hydrodynamics equations of a two-component fermionic superfiuid made of atoms of mass 
m are given by 
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—n + V ■ (nv) = (1) 
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^mv 2 + U(r) + fi(n, of) 
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where n(r, t) is the local density and v(r, t) is the local superfiuid velocity. Here n(r,t) = n^(r,t) + ni(r,t), with 
n-f(r, t) = ni(r,t) and v(r, t) = v-f(r, t) = v^(r, t). U(r) is the external potential and n(n, clf ) is the bulk chemical 
potential, i.e. the zero-temperature equation of state of the uniform system, which depends on the fermion-fermion 
scattering length a_p. The density n(r, t) is such that 

N = J n(r,t) d 3 r (3) 

is the total number of atoms in the fluid. In fact, due to the absence of the normal component, the superfiuid density 
coincides with the total density and the superfiuid current with the total current. 

Equations (fT]) and ||2} are nothing but the Euler equations of an inviscid and irrotational fluid. Since v is irrotational, 
it can be written as the gradient of a scalar field. The connection between superfiuid hydrodynamics and quantum 
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mechanics is made by the formula 



where 6{r, t) is the phase of the condensate wave-function 

5(r,t) = |S(r,t)| e»™ = (^ T (r,t)^(r,t)> , (5) 

with i/v(r, i) the fermionic field operator with spin component a =f, J, fl8l E9l|. Notice the factor 2m (Cooper pairs) 
instead of m in Eq. ® El El- 

The equations of superfluid hydrodynamics describe quite accurately static properties and low-energy collective 
modes of oscillation of wavelength A ^> £, where £ is the healing length of the superfluid. Recently, Combescot, Kagan 
and Stringari [2(| have suggested that 

(6) 

mVcr 

where v cr is the critical velocity of the Landau criterion for dissipation [H|2(j. According to Combescot, Kagan and 
Stringari [20| . in the BEC regime of bosonic dimers v cr coincides with the sound velocity, i.e. 



• (7) 

to on 



Instead, in the BCS regime v cr is related to the breaking of Cooper pairs through the formula 



^ = V m ' (8) 

where |A| is the energy gap of Cooper pairs Note that the equations of superfluid hydrodynamics {!]) and ^ do 
not take into account the effect of pair breaking. 

In the BCS-BEC crossover the bulk chemical potential of Eq. @ can be written as 

M(n, a F ) = £ (3^ 2 r7) 2/3 (/(„) - (9) 
where /(y) is a dimensionless universal function of the inverse interaction parameter 

y = ir— ( 10 ) 

where kp = (3n 2 n) 1 ^ 3 is the Fermi wavenumber (To| . One can parametrize /(y) as follows: 

f(y) = oti - OL 2 arctan ( a 3 y J , (11) 

where the values of the parameters aj, Qi2, ct3, Pi, 02, reported in Ref. [lOj], are fitting parameters based on asymptotics 
and fixed-node Monte-Carlo data |2l|. We will call in the following Monte-Carlo equation of state (MC EOS) the 
equation /i = fi(n, ap) obtained from ([9]) and (jlll) . 

Within the mean-field extended BCS theorv[19l. I22I ] , the bulk chemical potential fj, and the gap energy A of the 
uniform Fermi gas are instead found by solving the following extended BCS (EBCS) equations [22j, |23[ 



1 2 2m 1 / 2 Al/2 r , 2 / 1 
— = J A/ / d ^ - 2 - 



Al/2 / ^li-^Wr ( 12 ) 
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By solving these two EBCS equations one obtains the chemical potential /i as a function of n and ap in the full 
BCS-BEC crossover (see for instance Ref. [HI)- Note that EBCS theory does not predict the correct BEC limit: the 
molecules have scattering length a M — 2a F instead of a M — 0.6ap [Hj]. We call EBCS equation of state (EBCS EOS) 
the mean-field equation of state [i = fi(n,ap) obtained from Eqs. (|12p and (fT5|) . Of course, the MC EOS is much 
closer than the EBCS EOS to the MC results obtained in Ref. 
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III. SUPERFLUID NLSE FOR THE BCS-BEC CROSSOVER 



Inspired by the Ginzburg-Landau theory [241], by the density functional theory (DFT) [25|], and by the low-energy 
effective field theory (EFT) [H, we introduce the complex wave function 



*(r,t) 



n(r, t) 



which describes boson-like Cooper pairs with the normalization 

J |*M)| 2 d 3 r = ^ 



(14) 



(15) 



that is different from the normalization of the condensate wave function S(r, t) [23[ while the phase #(r, t) is the 
same [ill EH]- We now look for the simplest nonlinear Schrodinger equation of ^(r, t) which satisfies Eq. ([4} and 
reproduces the equations of superfluid hydrodynamics in the classical limit (U — *■ 0). We find [l,[ll[ that the nonlinear 
Schrodinger equation 



i*|*(r,*) 



- — V 2 + 2U(r) + 2/j(ra(r, t), a F ) 
4m 



*(r,t) 



(16) 



gives the equations of superfluid hydrodynamics in the classical limit (ft — ► 0). For finite ?L, this superfluid NLSE adds 
to the classical hydrodynamics equations a quantum pressure term 



T, 



QP 



8 m \Fn 



(17) 



containing explicitly the Planck constant fi (gradient correction in DFT, next-to-leading correction in low-energy 
EFT) 11]. Note that in the deep BEC regime from Eq. p^|) one recovers the familiar Gross-Pitaevskii equation for 
Bose-condensed dimers (molecules of two fermions) , where 



4-rrh 2 a dd {a F ) 
M(rw) = 2^ U 



(18) 



with addi^F) the dimer-dimer scattering length, which depends on the fermion-fermion scattering length a p. Within 
the mean-field theory one finds addi&F) — 2etp, whereas four-body theory and also MC many-body data yield 
add(ap) = 0.6a F [3. 



IV. DIRECT CURRENT JOSEPHSON EFFECT 



We use our time-dependent superfluid NLSE p6[) to study the direct-curent (DC) Josephson effect [HQ]. Consider 
a square-well barrier 



U(r) 



V for \z\ < d 
elsewhere 



which separates the superfluid into two regions, and assume a stationary solution 

*(r,t) - $(r) e 10 ^ e - l2flt/h 
with constant and uniform number supercurrent 

J = n(r)v(r) = 2$(r) 2 — V6»(r) . 

2m 

From the previous equation it follows (W) 2 = m 2 J 2 /(?i 2 <f >4 ) and also 



4m 



V 2 



m J 2 
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+ 2U(r) + 2fj,(n(T),ap) 



$(r) = 2/2 $(r). 



(19) 



(20) 
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(22) 
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FIG. 1: DC Josephson current J vs phase difference A8 for a superffuid Fermi gas at unitarity (y — 0), for three values of the 
energy barrier Vo- £f = h 2 (3n 2 n) 2 ^ 3 /(2m) is the Fermi energy. The width of the barrier is L — i/kp, where kp = (37r 2 n) 1//3 is 
the Fermi wave number. 



We solve this stationary superfluid NLSE by imposing a constant and uniform density n at infinity: 

*(r)-y| for |r|->co (23) 
Given <I>(r) at fixed J, the phase 0(r) is then obtained from 

e [r) = e M + f£ (24) 

The phase difference across the barrier is defined as 

A0 = 6(z = +oo) - 6{z = -co) . (25) 

This conditions allows us to establish the relationship between the current J and the phase difference A9. Figure 1 
reports our results obtained by using the bulk EBCS EOS at unitarity (ap = ±oo). As expected we recover the 
Josephson equation 

J = J sin(A0) (26) 

in the regime of high barrier (small tunneling, weak- link). In the limit of very small barrier (quasi-free transport, 
strong-link) Jo has its maximum value, given by 

Jo™" = n v cr (27) 

where v cr is the Landau critical velocity introduced in the previous section (20| . 

In Fig. 2 we plot J™ ax as a function of the inverse interaction parameter y = l/(kpaF), and compare our data 
(curves) with the results of Spuntarelli, Pieri and Strinati (symbols) [9] . Note that the data of Spuntarelli, Pieri and 
Strinati [9J have been obtained by solving the full self-consistent Bogoliubov-de Gennes equations of the EBCS theory 
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FIG. 2: Maximum Josephson current J™ ax vs inverse interaction parameter y = l/(kFa.F) in the BCS-BEC crossover. Solid 
curve: J™ ax based on pair breaking in the BCS regime [sL l20j] . Other curves: superfluid NLSE. Symbols: microscopic mean-field 
calculations of Ref. [9f. Four values of the energy barrier Vq/cf are considered: 0.025, 0.10, 0.2, 0.4. The width of the barrier 
is L = 4/k F - 



in the presence of the barrier flU]) . Figure 2 shows that the NLSE reproduces the DC Josephson results of Ref. 0], 
but only in the right side of the BCS-BEC crossover, i.e. from the unitarity limit to the deep BEC regime. It is not 
surprising that the NLSE fails in the BCS regime (y < 0), since it neglects completely the effect of pair breaking. The 
critical velocity v cr predicted by NLSE is the sound velocity c s also in the BCS regime. Nevertheless, it is remarkable 
that the NLSE works quite well in the BEC regime (y > 0) up to the unitarity limit (y = 0). 



V. ALTERNATE CURRENT JOSEPHSON EFFECT 



We now consider a high barrier (small tunneling, weak- link), without imposinga constant supercurrent J and discuss 
the alternate-current (AC) Josephson effect 0, 0] with superfluid Fermi atoms [8[. We start from the time-dependent 
superfluid NLSE (fll)|) and look for a time-dependent solution of the form 



*(r,t) = CA(t) $ A (r)+c B (t) $ B (r). 



(28) 



where $A(r) and $s(r) is the quasi-stationary solutions normalized to one and localized in region A and B respectively. 
In this way we obtainQ the following two-state model 



d 

ih—c A {t) = E A [t) c A (t) + K c B (t) 
d 

ih—c B (t) = E B {t) c B {t) + K c A (t) 



(29) 
(30) 



for the two complex coefficients c A (t) and c B (t), related to the number of atoms in the two regions. In our two-state 
model, E A (t) is the time-dependent energy in region A, given by 



E A (t) ~ / $ A (r) 



-— V 2 + 2U(r) + 2/i (2\c A (t)\H A (r) 2 , a F ) 



$ A (r) d 3 r. 



(31) 




FIG. 3: Zero-mode in the AC Josephson effect by solving AJJ equations. N — 1 6 40 K atoms between two symmetric regions of 
volume 25- 10 6 pm 3 , tunneling parameter K/ks = 10 -9 Kelvin and Fermi- Fermi scattering length = 1 /im, corresponding to 
y — 1.19. Solid line: population imbalance z(t); dashed line: phase difference <fi{t). Initial conditions: ip(0) = and 2(0) = 0.5 
(left); ip(0) = and z(0) = 0.999 (right). Adapted from Ref. [|. 



There is obviously a similar expression for the time-dependent energy E B (t). The constant coupling energy K describes 
instead the tunneling between the two regions A and B: 



K 



<5>a{y) 



-^V 2 + 2t/(r) 



$s(r)d 3 r. 



(32) 



From our previous analysis of the DC Josephson effect, we expect that this expression is correct only in the right 
side of the BCS-BEC crossover. However, to extend the study the Josephson effect to the left side of the BCS-BEC 
crossover one may use K in Eqs. ([29]) and (|30[) as a phenomenological parameter a la Feynman pH, [l7j. 
We can write the complex coefficient ca (t) as 



c A {t) 



N A (t) 



iS A (t) 



(33) 



with NA(i) number of atoms and 9(t) phase in region A. Again, a similar expression holds for cb(£). In terms of the 
phase difference 



and relative number imbalance 



z(t) 



cp(t) = e B {t)-e A (t) 

N B (t)-N A (t) _ N B (t)-N A (t) 
N A (t)+N B {t) N 



the two-mode equations give 



2K 



z{t) = -—^Jl-ziif sin ¥>(*). 



<p(t) 



N 



2V 



N 



2V 



2K 



z{t) 



cos tp(t) . 



(34) 
(35) 

(36) 
(37) 



These are the atomic Josephson junction (AJJ) equations describing the oscillations of N Fermi atoms tunneling in 
the superfluid state between region A and region B, of equal volume V Q. Notice that equations (|36[) and l|37p 
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FIG. 4: 7r-mode in the AC Josephson effect by solving AJJ equations. All physical parameters are the same as in Fig. 3. 
Upper panel: population imbalance z(t); lower panel: phase difference ip{t). Initial conditions: tp(0) = and z(0) = 0.5 (solid); 
ip(0) = and z(0) = 0.999 (dashed). 



generalize the BJJ equations obtained by Smerzi et al. [3J] for Bose-Einstein condensates. Moreover, these equations, 
linking the tunneling current 

N KN 



1 V IY I V / / 

J = —z— = — v 1 — z 2 sin tp — Jo v 1 — z 2 sin tp (38) 

to the phase difference tp, reduce to the familiar Josephson expression J = J$ sia(tp) in the appropriate limit \z\ <C 1. 
The nonlinear AJJ equations can be linearized around the stable stationary solution 

z = and tp — 2nj (39) 

where j is an integer. In this way one finds the following frequency of small oscillation 



»o = ^ + ^f (40) 

which is called zero-mode. Here c s is the sound velocity computed at the mean density n = N/V of the supcrfiuid. 
This zero-mode is the analog of the Josephson plasma oscillation in superconducting junctions 0. Figure 3 reports 
the zero-mode oscillations of z(t) and tp(t), in the case of 40 K atoms setting aj? = 1 /im. The oscillation starting 
from z(0) = 0.5 indicates that the solution flU)) of the linearized equations (|31)|) and (|3"?| are fairly accurate even for 
finite and not quite small amplitude. Eventually however, for very large amplitude, z(0) = 0.999, deviations from the 
harmonic approximation become quite visible. 

We observe that the AJJ equations (|36|) and 1|37|) produce also a 7r-mode solution with z = and tp = n(2j + 1) and 
the self-trapping solution with population imbalance (z ^ 0) (for details see [8j). Figure 4 reports 7r-mode oscillations 
for two different initial unbalance. Note the significant frequency increase and waveshape distortion induced by 
nonlinear effects. 



VI. CONCLUSIONS 



This paper reviews a superfluid NLSE providing the hydrodynamic equations of Fermi superfluids plus a gradient 
correction. Both hydrodynamics equations and superfluid NLSE are reliable to investigate static properties and low- 
energy collective dynamics. It is important to stress that the equations of superfluid hydrodynamics are nothing else 
than the time-dependent local density approximation (LDA) for an irrotational system. The advantage of using the 
NLSE is that can take into account surface and shape effects beyond LDA, and these can be relevant for a small 
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number of particles [2Sj . We have shown that in the study of the DC Josephson effect, the superfluid NLSE works 
quite well at the right side (BEC regime) of the BCS-BEC crossover up to the unitarity limit of infinite scattering 
length. In addition, we have suggested that for studying DC and AC Josephson effects in weakly-linked atomic 
superfluids the NLSE can be used in the full BCS-BEC crossover within the two-mode approximation. In this case 
the tunneling energy coefficient must be taken as a phenomenological parameter at the left side (BCS regime) of the 
crossover. 

We thank Andrea Spuntarelli and Pierbiagio Pieri for making available their data @. LS acknowledges a reseach 
grant from GNFM-INdAM. This work has been partially supported by Fondazione CARIPARO. 
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